Stability for Walls in Ferromagnetic Nanowire 



O ■ G. Carbou 1 and S. Labbe 2 

o i 

s ! ' 1 MAB, Universite Bordeaux 1, 351 cours de la Liberation, 33405 Talence cedex 

carbouOmath . u-bordeauxl . f r 
2 Universite Paris Sud, Orsay labbe@math.u-psud.fr 



o 

(N 

< 

c3 



Summary. We study the stability of travelling wall profiles for a one dimensional 
model of ferromagnetic nanowire submitted to an exterior magnetic field. We prove 
that these profiles are asymptotically stable modulo a translation-rotation for small 
applied magnetic fields. 



1 Model for ferromagnetic nanowires 

Ferromagnetic materials are characterized by a spontaneous magnetization 
described by the magnetic moment u which is a unitary vector field linking 
the magnetic induction B with the magnetic field H by the relation B = H+u. 
The variations of u are described by the Landau-Lifschitz Equation 
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— = -u A H e - u A (u A H e ) (1) 

VO ' 

. where the effective field is given by H e = Au + hd{u) + H a , and the demag- 

""p^ ' netizing field hd{u) is deduced from u solving the magnetostatic equations: 

> 

div B = div (H + u) = and curl H = 

where H a is an appplied magnetic field. 

For more details on the ferromagnetism model, see [2], [9], [14] and [18]. 
For existence results about Landau Lifschitz equations see [3], [4], [10], [17]. 
?—i ' For numerical studies see [8], [12] and [13]. For asymptotic studies see [1], [5], 

[7], [15] and [16]. 

In this paper we consider an asymptotic one dimensional model of ferro- 
magnetic nanowire submited to an applied field along the axis of the wire. We 
denote by (ei, 62,63) the canonical basis of IR 3 . The ferromagnetic nanowire 
is assimilated to the axis lRe\. The demagnetizing energy is approximated by 
the formula hd(u) = —U2&2 — 1*363 where u = (ui, 112, U3) (this approximation 
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of the demagnetizing energy for a ferromagnetic wire is obtained using a BKW 
method by D. Sanchez, taking the limit when the diameter of the wire tends 
to zero in [16]). We assume in addition that an exterior magnetic field <5e x is 
applied along the wire axis. 

To sum up we study the following system 

— = — u A h$(u) — u A (u A h$(u)) 
d 2 u 

with h s (u) = - u 2 e 2 - u 3 e 3 + Sei 

For 5 = 0, that is without applied field, we observe in physical experi- 
ments the formation of a wall breaking down the domain in two parts: one 
in which the magnetization is almost equal to e\ and another in which the 
magnetization is almost equal to —e\. Such a distribution is described in our 
one dimensional model by the following profile M : 



Mo = . (3) 




This profile is a steady state solution of Equation (2) with 6 = 0. We 
prove in [6] the stability of the profile M for Equation (2) without applied 
field (when 5 = 0). 

When we apply a magnetic field in the direction +ei (that is with S > 0) 
since Landau-Lifschitz Equation tends to align the magnetic moment with 
the effective field, we observe a translation of the wall in the direction — e\. 
Furthermore, we observe a rotation of the magnetic moment around the wire 
axis. This phenomenon is described by the solution of (2) 

U 5 {t,x)=R St (M (x + St)) (4) 
where Rg is the rotation of angle 8 around the axis Me\\ 

(10 \ 
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We study in this paper the stability of Us, we prove that for a small 6, Us 
is stable for the H 2 norm and asympotically stable for the H 1 norm, modulo a 
translation in the variable x and a rotation around JRe\. This result is claimed 
in the following theorem: 

Theorem 1. There exists <5o > such that for all 5 with |<5| < So then for 
e > there exists r\ > such that if \\u(t = 0,x) — Us(t = 0, ir)|| h 2 < V then 
the solution u of Equation (2) with initial data u(t = 0,x) satisfies: 
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Vt > 0, \\u(t,x) - U s {t,x)\\ H 2 < e. 

In addition there exists and 9^ such that 

\\u(t,x) - R gao {Us{t,x + £7oo))||i?i — ► when t — ► +oo. 

This result is a generalization of the stability result concerning the static 
walls when S = in [6]. It looks like the theorems of stability concerning 
the travelling waves solutions for semilinear equations like Ginzburg Landau 
Equation (see Kapitula [11]). Here we have three new difficulties. The first 
one is that the magnetic moment takes its values in the sphere and not in a 
linear space. In order to work with maps with values in a linear space we will 
use a mobile frame adapted to Landau-Lifschitz equation and we will describe 
in Section 2 the magnetic moment in this mobile frame. The second difficulty 
is that we have here a two dimensional invariance family for Equation (2) 
whereas Ginzburg Landau-Equation is only invariant by translation. This is 
the reason why we must use in the perturbations description the translations 
and the rotations (see Section 3). The last difficulty is that Landau-Lifschitz 
Equation is quasilinear, and then we have to couple variational estimates and 
semi-group estimates to control the perturbations of our profiles. Section 4 is 
devoted to these estimates. 



2 Landau-Lifschitz Equation in the mobile frame 
2.1 First reduction of the problem 



For u a solution of Landau-Lifschitz Equation (2) we define v by v(t, x) = 
Rst(u(t,x — St)) (that is u(t, x) = Rst(v(t,x + St))). A straightforward cal- 
culation gives that u satisfies (2) if and only if v satisfies 



r dv 

di 



= -v A h(v) - v A (v A h(v)) 

d 2 v 



(5) 



h(v) = — - v 2 e 2 - v 3 e 3 



In addition Us is stable for (2) if and only if M is stable for (5), that is 
we are led to study the stability of a static profile, which is more convenient. 



2.2 Mobile frame 



Let us introduce the mobile frame (Mq(x), Mi(x), M 2 ), where 



M x (x) 



{ J_ 

chx 


\-th.T, 



(T 

and M-2 --- | 1 
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Let v : 2R+ x JR X — ► S 2 C t 3 be a little perturbation of M . We can 
decompose v in the mobile frame writting 

v(t,x) = r 1 (t,x)M 1 (x)+r 2 {t,x)M 2 + yjl - r\ - r|M (x). 

Now we can obtain a new version of Landau-Lifschitz Equation: v satisfies 
(5) if and only if r = (n, r 2 ) satisfies 

| = (£ + «)r + G(r)(g) + H(x,r,^) (6) 



where 

I ( 1 ■. I 1 ■ lllM'VMllll' / 1 o l \ ■( il l 1 i\ ■ /' / / u 'l 1 1 l J | 



• the linear operator C is given by C — JL with J = I ^ j I and 



d 2 

L = + 2th 2 x- f 



• the linear perturbation due to the presence of the applied magnetic field 

d 

Sei is given by 51 with I = — — h thx, 

ox 

• the higher degree non linear part is G{r){^ 1 ), where G{r) is a matrix 
depending on r with G(0) = 0, 

• the last non linear term H (x, r, is at least quadratic in the variable 

dv 

In addition the stability of the profile M for Equation (5) is equivalent to 
the stability of the zero solution for Equation (6). 



3 A new system of coordinates 

We remark that L is a self adjoint operator on L 2 (M), with domain H 2 (M). 

d 

Furthermore, L is positive since we can write L = I* o I with I = — — h thx, 

ox 

and Kcr L is the one dimensional space generated by — — . 

chx 

The matrix J being invertible, Ker C is the two dimensional space gener- 
ated by v\ and v 2 with 

( \ 

v\(x) = _\_ I , v 2 (x) = 

V chx ' 

We introduce £ = (Ker We denote by Q the orthogonal projection 
onto £ for the L 2 (1R) scalar product. 
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Landau-Lifschitz equation (5) is invariant by translation in the variable x 
and by rotation around the axis e\. Therefore for A = (6, a) fixed in M 2 , M A 
defined by M A (x) = Rg(M (x—a)) is a solution of Equation (5). We introduce 
Ra(x) the coordinates of Ma(x) in the mobile frame (Mi(x), M 2 (x)): 



R A {x) 

The map \P given by 



\M A {x)-M 2 



<P : M 2 x £ — > H 2 (M) 

(A, W) i — ► r(x) = R A (x) + W(x) 

is a diffeomorphism in a neighborhood of zero. Thus we can write the solution 
r of Equation (6) on the form : 

r(t,x) = R A(t) {x) + W(t,x) 

where for all t, W(t) € £ and where A : R+ ^ R 2 . 

We will re- write Equation (5) in the coordinates (A, W). Taking the scalar 
product of (5) with v\ and v 2 we obtain the equation satisfied by A, and using 
Q the orthogonal projection onto £, we deduce the equation satisfied by W. 
After this calculation we obtain that r is solution of Equation (5) if and only 
if (A, W) satisfies the following system 



dW d 2 W dW 

— = (£ + 51 + K A )W + H^x, A, W)(-^-) + K 2 (x, A, W, —) 

dA KAnxr dW „ 
Tt =M{W,—,A) 



(7) 



where 

• K A : H 2 (M) — > £ is a linear map satisfying 

3K U VAelR 2 ,VWe£, \\K A W\\ L , {m < ^i|yl|||VK||^( K ) (8) 

• the non linear terms take their values in £ and satisfy that there exists a 
constant K 2 such that for \A\ < 1 and for all W £ £ 

d 2 W 

||ftiM, W)(-^-)\\ L 2 {m < K 2 \\W\\ Hl{m \\W\\ H 2 {m 

(9) 

\\n 2 (.,A,w,^)\\ H1{R) <K 2 \\w\\ 2 HHm 

• M:H 1 (M)x L 2 (M) x M 2 — ► M 2 satisfies 

dW 

3K 3 , V/l such that \A\ <l,VWe£, \M(W, A)\ < K 3 \\W\\ H i {m 

(10) 
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Theorem 1 is equivalent to the following Proposition: 

Proposition 1. There exists <5 > such that for 5 with \S\ < So, we have the 
following stability result for Equation (7): for e > there exists r\ > such 
that if\A \ < f] and i/||Wo||# 2 < f] then the solution (A, W) of (7) with initial 
value (A ,W ) satisfies 

1. for all t > 0, \\W{t)\\ H 2 < e and \A\ < e, 

2. \\W{t)\\}ji tends to zero when t tends to +oo, 

3. there exists A^ e M 2 such that A(t) tends to A x when t tends to +oo. 
The last section is devoted to the proof of Proposition 1 . 

4 Estimates for the perturbations 
4.1 Linear semi group estimates 

On £ we have Re (sp C) c] — oo, — 1]. In particular this fact implies that the 
H 2 norm is equivalent on £ to the norm \\Cu \\ L 2. Furthermore it implies good 
decreasing properties for the semigroup generated by C. We first prove that 
this decreasing property is preserved for the linear part of the Equation on W 
in (7) for a little applied field, and if we assume that A remains little. 

The operator I is an order one operator dominated on £ by £, thus there 
exists <5o > such that if \8\ < So, Re (sp C + SI) c] — oo, — 1/2[. 

Let us fix S such that |<5| < <5o- With Estimate (10), if A remains small, JCa 
is a little perturbation of C + SI. This implies that for A little, the semigroup 
generated by C + SI + QJCa has the same good decreasing properties than £, 
that is there exists v$ > such that if \A(t)\ remains less than vq for all t, 
then there exists K4 and [3 > such that 



We can then use the Duhamel formula to solve the equation on W in (7): 



and then using the estimates (9) and (11) we obtain that if |^l(i)| remains less 
than vq then there exists K 5 such that 

f t e -0(t-s) 

^ -flfllTFT II / C II,,,/ Nil II,,,/ 111 



\\S A (t)W \\ H i < K 4 e~^ || Wo \\ H i 
P -0t 

<K 4 __||Wb|| L 2. 



(11) 





(12) 




Stability for Walls in Ferromagnetic Nanowire 



7 



4.2 Variational estimates 

We see that Estimate (12) is not sufficient to conclude since we have the H 2 
norm of W in the right hand side of this estimate. In order to dominate this 
H 2 norm, we multiply the equation on W in (7) by J 2 C 2 W and we obtain 
that there exists a constant K & : 

j t \\LW\\ 2 L2 + \\L-*W\\\ 2 (1 - K 6 \\LW\\ L2 ) < 

; From this estimate we deduce that if ||LVF|| L 2 < ±,th.enl-K 6 \\LW\\ L 2 
is positive, thus ^H-^W^I^ is negative and ||LVF|| L 2 remains less than 
So if \\LW \\ L 2 < -U, then for all t \\LW(t)\\ L 2 < \\LW \\ L *. This property 
gives a bound for the H 2 norm of W since the H 2 norm is equivalent on £ 
to ||LW||i,2, and reducing the H 2 norm of Wo, we obtain the first part of the 
conclusion 1 in Proposition 1. 



4.3 Conclusion 

Let us assume that ||I<Wo||i,2(jj) < Then for all t, \\W(t)\\ h 2 (m) < 

Ci||iVF(t)|| L 2 < ||-LWo||z,2 < C 2 ||Wo||jj2 ( b), where d and G 2 are constants. 
Multiplying (12) by (1 + t) 2 , defining G{t) = max(l + s) 2 \\W{s)\\ H i , we 

obtain that there exists a constant K 7 such that if \A(t)\ remains less than vo 
we have: 

G(t) < K 7 G(0) + K 7 G(t)\\W \\ H 2 + K 7 (G(t)) 2 
If we suppose in addition that ||Wo||jf2 < — we obtain that 

< K 7 G(0) - ^G(i) + K r (G(t)) 2 := P(G(t)) (13) 

The polynomial map P(£) = K 7 £ 2 - ±£ + K 7 G(0) has for G(0) small 
enough two positive roots. We denote by £(G(0)) the smallest one. For G(0) 
little enough we have G(0) < £(G(0)) < 2K 7 G(0) (we can a priori assume 
that K 7 > 1 for example). Estimate (13) implies that for all t, G(t) < £(G(0)) 
that is 

Vt > 0,\\W(t)\\ m{R) < — — < —. (14) 

This implies that || W^WIIff^B) tends to zero when t tends to +00. It remains 
to prove that A remains less that vq and admits a limit when t tends to +00. 

Plugging Estimate (14) in the equation on A in (7) and using (10), we 
obtain that ^ is integrable on M + , that is A admits a limit when t tends to 
+00. Furthermore, by integration we have 

V t, \A(t)\ < \A(0)\ + f K 3 2I ^ 7G{ ^ ds < \A(0)\ + ttK 3 K 7 G(0) 
Jo 1 + s 

Reducing \A \ and G(0) = || WoWmiR) wc obtain that for all t, \A(t)\ remains 
less than i/q, which justifies all our estimates a posteriori. 
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